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ABSTRACT 

The abelian sandpile models feature a finite abelian group G generated by the 
operators corresponding to particle addition at various sites. We study the canonical 
decomposition of G as a product of cyclic groups G = x Z^,^ x ■ ■ ■ x Zfi^ where g 
is the least number of generators of G, and di is a multiple of di+i. The structure of G 
is determined in terms of the toppling matrix A. We construct scalar functions, linear 
in height variables of the pile, that are invariant under toppling at any site. These 
invariants provide convenient coordinates to label the recurrent configurations of the 
sandpile. For an L x L square lattice, we show that g = L. In this case, we observe 
that the system has nontrivial symmetries, transcending the obvious symmetries of 
the square, viz. those coming from the action of the cyclotomic Galois group Gal/, of 
the 2(L + l)-th roots of unity (which operates on the set of eigenvalues of A). These 
eigenvalues are algebraic integers, whose product is the order With the help of 
Gall, we are able to group the eigenvalues into certain subsets whose products are 
separately integers, and thus obtain an explicit factorization of \G\. We also use Gal^ 
to define other simpler, though under-complete, sets of toppling invariants. 
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1. Introduction 

The concept of self-organized criticality was proposed by Bak, Tang and Wiesenfeld 
in 1987, and has given rise to growing interest in the study of self-organizing systems. 
Bak et al argued that in many natural phenomena, the dissipative dynamics of the 
system is such that it drives the system to a critical state, thereby leading to ubiquitous 
power law behaviors [1,2]. This mechanism has been invoked to understand the power- 
laws observed in turbulent fluids, earthquakes, distribution of visible matter in the 
universe, solar flares and surface roughening of growing interfaces [3-5] . 

Sandpile automata are among the simplest theoretical models which show self- 
organized criticahty. A specially nice subclass consists of the so-called abehan sandpile 
models (ASM's) [6]. There have been many numerical as well as analytical studies of 
the ASM. The case when there is a preferred direction of particle transfer turns out 
to be equivalent to the voter model, and all the critical exponents can be determined 
in all dimensions [7]. When there is no preferred direction, the model turns out to be 
related to the g — > limiting case of the Potts model [8]. The problem has been solved 
exactly in the mean field limit [9-11]. In two and three dimensions, only some of the 
critical exponents of the problem are known analytically [6,8,12]. 

Most of these papers have been concerned with the critical properties of the ASM 
in the thermodynamic limit of large system sizes. In this paper, our interest is rather 
in the properties of the finite automata. There have been only a handful of papers 
addressing these so far (and mainly in two dimensions). Creutz' paper [13] exhibits 
very interesting examples of geometrical patterns displayed by a special configuration 
(the so-called identity configuration) of the ASM on a square lattice. Liu et al [14] 
have studied the patterns obtained by relaxing some simple unstable configurations. 
Wiesenfeld et al [15] have studied the periods of deterministic ASM, again on a square 
lattice. These studies have been extended by Markosova and Markos [16] to lattices 
of size up to 19. This paper was inspired partly by the 'experimental' results of these 
authors. We will concentrate on properties such as the structures of the abelian group 
and of the space of recurrent configurations. 

The plan of this paper is as follows. In section 2, we give the definition of general 
ASM, and review their basic properties. In section 3, we construct, for a general ASM, 
a set of functions, defined on the space of recurrent configurations, which are invariant 



under topplings, and which can be used to label those configurations. With these 
functions, we show how to determine the structure of the abelian group of the ASM 
in terms of the normal form decomposition of its toppling matrix. In section 4, we 
consider the problem of computing the rank of the abelian group when the ASM is 
defined on an Li x L2 rectangular portion of the square lattice, and show that the rank 
equals L on an L x L square. In section 5, we recapitulate briefiy the results from 
general Galois theory needed by us, and use them to study the Galois group of the 
characteristic polynomial of the toppling matrix, and construct, in section 6, another 
set of algebraic functions, invariant under topplings, thereby extracting information 
about the abelian group from Galois theory. This information is incomplete and, we 
show that the method does not give the full structure of the group in general. In 
section 7, we remark on two other properties of the ASM on a square lattice. Firstly, 
we derive a sharper upper bound on the time period of the deterministic ASM studied 
in [15]. The second one is an interesting connection between the structure of identity 
configurations on 2L x 2L and (2L + 1) x (2L + 1) lattices. Some technical material 
involving detailed calculations is given in appendices A and B, while in appendix G we 
display some identity configurations for some small square lattices. In appendix D, we 
show how the toppling invariants discussed by Lee et al [17] are a special case of the 
invariants discussed in this paper. 



2. Preliminaries and notation 

The general ASM is defined on a set of N sites, labelled say by integers 1 to A^. 
Each site i is assigned an integer variable Zi, called the height of the sandpile at site i. 
The time evolution of the sandpile is defined in terms of the following two processes: 

(1) Addition of particles: We choose a site at random, and increase its height by 1, 
while the heights at other sites remain unchanged. The probability of choosing 
the k^^ site is denoted by p^. For simpficity we take all p^s nonzero. [This 
condition is needed to ensure the uniqueness of the steady state.] 

(2) Relaxation: If the height at some site j equals or exceeds a prespecified value 
ZcritU), it is said to be unstable and topples, loosing some grains of sand which 



either fall on other 'neighbouring' sites (whose height increases as a consequence), 
or drop out of the system. The updating of heights is specified in terms of an 
N X N integer matrix A, called the toppling matrix and satisfying A^, > and 
Aij < for i 7^ j. If there is a toppling at some site j, all heights {zi} are 
updated according to the rule 

If Zj > Zcritij), then Zi ^ Zi — Ay for all i. (2.1) 

We may assume that 

Zcritij) = ^jj: (2-2) 

so that the values of Zj in a stable configuration are between and Ajj — 1. 

A toppling at one site may make other sites unstable. A sequence of topplings 
caused by adding a single particle is called an avalanche. When all the unstable sites 
have toppled, we are left with a new stable configuration. This defines a single step of 
time evolution. At the next time, a new particle is added at a random site, and the 
system is allowed to relax, and so on. 

It is convenient to define operators ak, k = 1 to N, corresponding to the process of 
adding a particle at site k, and allowing the system to relax. These operators map 
stable configurations into stable configurations. 

A general analysis of abelian sandpile models was carried out in [6]. It was 
shown that the operators {aj.} commute with each others. This allows a simple 
characterization of the critical steady state: Only a small subset of all stable 
configurations occur with non-zero probability in the steady state. These are called 
recurrent configurations. Their number is equal to Det A, and in the steady state they 
occur with equal probability. The a^'s map the space TZ of recurrent configurations 
onto itself, and are invertible on TZ. 

Let G be the abelian group generated by the operators {a^}, i — 1 to N. This is a 
finite group as these operators satisfy the closure relations [6] 

N 

Haf'^^I, for all J. (2.3) 

1=1 

The order of G, denoted by is equal to the number of recurrent configurations. 
This is a consequence of the fact that if C and C' are any two recurrent configurations. 



then there is an element g & G such that C = gC . We thus have 

IGj = |7e| = DetA. (2.4) 



3. Toppling invariants and the group structure for general ASM 

The space of all configurations {zi} (with non-negative heights Zi) constitutes a 
commutative semigroup over the given set of N sites, with the operation given by 
sitewise addition of heights followed by relaxation if necessary. One can define an 
equivalence relation on this semigroup by declaring two configurations {zi} and {z[} 
equivalent (under toppling) if and only if there exist N integers rij such that 

z[ — Zi — ^ij'i^i 1 for alH. (3.1) 

3 

Each equivalence class contains one and only one recurrent configuration. Indeed to 
any configuration {^i}, one can associate a recurrent configuration C by letting 

C[zi]^Y[afC\ (3.2) 

i 

where C* is any fixed recurrent configuration. Then if {zi} and {z[} are related as in 
(3.1), one easily checks that C[zi\ = C\z'^ on account of the relations (2.3). (Note that 
with the definition (3.1), two stable configurations may be equivalent under toppling.) 

Toppling invariants are scalar functions defined on the space of all configurations of 
the sandpile, such that they take the same value for any two configurations which are 
equivalent under toppling. Toppling invariants which are linear in the height variables, 
and which are conserved modulo various integers were first introduced by Lee and 
Tzeng [17] in the context of specific deterministic one-dimensional ASM. Their results 
will be shown to be particular cases of our general construction (see Appendix D). 



Given the toppling matrix A, we define N rational functions Qi {i = 1 to N) by 
setting 

Qi({zj})^Yl^^'j (3.3) 
j 

It is straightforward to prove that the functions Qi are toppling invariants. Indeed, 
under toppling at site k, the configuration C = {zj} changes to C = {z'j — zj — Aj/j}, 
and from the linearity of the functions Qi in the height variables, we have 



Qi{C') = Qi{C) - J2 ^^'^jk = Qi{C) mod 1. 



(3.4) 



The functions Qi take rational values, but they are easily made integer-valued upon 
multiplication by some adequate integer. Being toppling invariants, the functions Qi 
can be used to label the recurrent configurations, and thus the space TZ can be replaced 
by the set of A'^-uples (Qi, Q2, • • • , Qn)- However, the labelling by the QiS is generally 
overcomplete, they not being all independent. A simple example will readily establish 
that. 

Throughout this paper, we shall consider in detail the special case of the ASM 
defined on an Li x L2 rectangular portion of a two-dimensional square lattice. We 
choose the toppling matrix to be the discrete Laplacian, whose diagonal entries are 
given by An — 4, and the off-diagonal entries Ay = — 1 or according to whether 
the sites i and j are nearest neighbours or not. Note that this implies open boundary 
conditions on all four boundaries of the rectangle, and thus any toppling there involves 
a loss of sand. 

To be specific, consider the case L1 — L2 — 2, with the configurations specified as 

Zl Z2 



In this case A is a 4 x 4 matrix of determinant 192, and we find 
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The definition (3.3) yields the four invariants 

1 



(3.5) 



24 



{7zi + 2z2 + 2z2, + z^) mod 1, 



(3.6a) 



^2 = ^(2^1 + 7^2 + -^3 + 2^4) modi, (3.66) 

Q3 = ^{2zi + Z2 + 7z3 + 2z4) modi, (3.6c) 

Q4 = ^ (-21 + 2^2 + 2^3 + 7^4) modi. {3.6d) 

They satisfy three hnear relations 

Qi + Q2 = mod 1/8, (3.7a) 
Q3 = ^Qi-Q2 modi, (3.76) 
Q4 = -Qi + 4Q2 modi, (3.7c) 

from which one sees that onfy two independent invariants remain, which we choose 
integer-valued for convenience 

h = 24Qi = 7zi + 2^2 + 2zs + Z4 mod 24, (3.8a) 
h = S{Qi + Q2) = 3zi + 3^2 + Z3 + Z4 mod 8. (3.86) 

These two invariants provide a complete labelling for the space of recurrent 
configurations, of cardinality 192. 

As this very simple example already shows, for an arbitrary N x N matrix A, we 
construct N toppling invariants Qi, but they are generally not independent. It thus 
seems desirable to isolate a minimal set of independent invariants, as we did above to 
obtain the mininal set given by (7i,/2)- We now show how this can be done for an 
arbitrary ASM using the classical theory of Smith normal form for integer matrices. 

In current mathematics literature, as in Jacobson [18], this is often discussed for 
matrices with entries in a principal ideal domain (the ring of integers being a prime 
example). Prom Theorems 3.8 and 3.9 of Jacobson, any nonsingular N x N matrix A 
can be written in the form 

A = ADB, (3.9) 

where A and B are N x N integer matrices with determinants ±1, and D is a, diagonal 
matrix 

Dij = didij (3.10) 

such that 



1. (ij is a multiple of (ii+i for alH = 1 to — 1, and 

2. di — Ci-i/ei , where stands for the greatest common divisor of the determinants 
of all the {N — i) X {N — i) submatrices of A (set ejv = 1). 

Thus the matrix D is uniquely determined by A, but the matrices A and B are far 
from unique. The integers d-i are are called the elementary divisors of A. An efficient 
algorithm to compute them can be found in Cohen [19]. 

In terms of the decomposition (3.9), we define the set of scalar functions Ii{C) by 



(3.11) 



Due to the unimodularity of A, these functions are integer- valued. As argued for the 
Q's in Eq. (3.4), we see using (3.9) that Ii{C) so defined are invariant under the 
toppling of any site. Clearly, only those invariants li with di^ 1 are nontrivial. Note 
that they are written in terms of A, and hence not unique. It will be obvious from the 
discussion below that the set of non-trivial U is always minimal and complete. 

As suggested by the notation used in the above example, this second set {U} 
is a minimal set of invariants chosen from the overcomplete set {Qi\. The Smith 
decomposition precisely shows how to combine the overcomplete invariants Qi so as to 
obtain a complete set. Indeed, it is easily checked that the /j's can be written in terms 
of the QiS as 

li — di BijQj mod di. (3.12) 



For instance, in the above 2x2 example, one finds the following Smith decomposition 
(written in the form A~^^ — DB) 



/ 7 2 2 1\ 

3 3 11 

-10 

V-2 -1 0/ 



diag(24,8,l,l) 



/ 1 0\ 
110 
-4110 
V-7 2 -2 1/ 



(3.13) 



From the first two rows of A ^ and (3.11), one recovers the two invariants Ji, I2 of 
(3.8). From the matrix B and (3.12), one finds the relations (3.8) between {h^h) 



and (Qi, Q2), while the trivial character of Is, I4 leads to the linear relations (3.7b-c) 
among the Q's. 

Let us now show that the set not only forms a complete set of toppling 
invariants, but that they also determine the structure of the abelian group G. 

Let g be the number of rf^ > 1. With each recurrent configuration of the ASM, we 
associate a g-uple {Ii, I2, ■■■■Ig), where < /j < di. The total number of distinct g- 
uples is nf=i — 1^1 (since A and B in the Smith decomposition of A are unimodular). 

We first show that this mapping from the set of recurrent configurations to g^-uples 
is one-to-one. Let us define operators by the equation 

N 

ei = l[af'\ l<i<g. (3.14) 

Acting on a fixed recurrent configuration C* = {zj}, 6% yields a new recurrent 
configuration, equivalent under toppling to the configuration {zj + Aji). If the g- 
uple corresponding to C* is [I^, /|, . . . , /*), it is easy to see from (3.11) that CjC* gives 
a configuration whose toppling invariants are //^ = + ^j^. By operating with these 
operators {cj} sufficiently many times on C*, all \G\ values for the s^-uple (Ji, I2, . . . , Ig) 
are obtainable. Thus, there is at least one recurrent configuration corresponding to any 
51-uple {Ii, I2, . . . , Ig). As the total number of recurrent configurations exactly equals 
the total number of g^-uples, we see that there is a one-to-one correspondence between 
the 51-uples (/i, h: - ■ ■ 1 Ig) and the recurrent configurations of the ASM. 

To express the operators aj in terms of e^, we need to invert the transformation 
(3.14). This is easily seen to be 

9 _^ 

aj = Yl , 1 < J < iV. (3.15) 

1=1 

Thus the operators generate the whole of G. Since Cj acting on a configuration 
increases /j by 1, leaving the other invariants unchanged, and since /j is only defined 
modulo di, we see that 

ei"^' = /, for i = 1 to g. (3.16) 

Note that the definition (3.14) for Cj makes sense for i between g+1 and N, and implies 
relations among the aj operators (cj = / from (3.16)). 



This shows that G has the canonical decomposition as a product of cychc groups 



G = Zd,x Zd,x ...X Zd, 



(3.17) 



with the diS defined in (3.10). We thus have shown that the generators and the group 
structure of G for an arbitrary ASM can be entirely determined from its toppling 
matrix A, through its normal from decomposition (3.9). 

The invariants {U} also provide a simple additive representation of the group G. 
As mentioned at the beginning of this section, one can define a binary operation of 
"addition" (denoted by ®) on the space TZ of recurrent configurations by adding heights 
sitewise, and then allowing the resultant configuration to relax. Prom the linearity of 
the liS in the height variables, and their invariance under toppling, it is clear that 
under this addition of configurations, the U also simply add, i.e. for any recurrent 
configurations Ci and C2, one has 



The Jj's provide a complete labelling of TZ. There is a unique recurrent configuration, 
denoted by Cj^, for which all Ii{C\^) are zero. Also, each recurrent C has a unique 
inverse — C, also recurrent, and determined by Ii{—C) — —Ii{C) mod dj. Therefore 
the addition ® is a group law on TZ, with identity given by C\^. A simple recursive 
algorithm to compute has been given by Creutz [13]. 

There exists a one-to-one correspondence between the recurrent configurations of 
ASM, and the elements of the group G: we associate with the group element g & G, 
the recurrent configuration gC^^. It is then easy to see from (3.18) that for all g,g' & G 



Thus the recurrent configurations with the operation © form a group which is 
isomorphic to the multiplicative group G, a result first proved in [13]. 

The invariants {li} provide a simple labelling of all the recurrent configurations. 
Since a recurrent configuration can also be uniquely specified by the height variables 



Ii{Ci®C2) = Ii{Ci) + Ii{C2) mod di. 



(3.18) 



gCid © g'Cid = {gg')Cid . 



(3.19) 



{zi}, the existence of forbidden subconfigurations in ASM's implies that these heights 
satisfy many inequaUty constraints. [Certain patterns cannot appear inside any 
recurrent configuration, and hence are called forbidden subconfigurations [6]; for 
instance, two neighbouring sites can never both have zero heights.] 

4. The rank of G for a rectangular lattice 

For a general matrix A, it is difficult to say much more about the group structure 
of G. In the rest of this paper, unless otherwise specified, we shall consider the special 
case when A is the toppling matrix corresponding to a finite Li x L2 rectangle of the 
(two-dimensional) square lattice. Here it is more convenient to label the sites not by 
a single index i going from 1 to = L1L2, but by two Cartesian coordinates {x,y), 
where 1 < x < Li and 1 < y < ^2- The toppling matrix is the discrete Laplacian, 
given by A{x, y; x, y) = 4, A{x, y; x', y') = —1 if the sites (x, y) and {x', y') are nearest 
neighbours {i.e. \x — x'\ + \y — y'\ = 1), and zero otherwise. Without loss of generality, 
we can assume that Li > L2. 

The relations satisfied by the particle addition operators a{x, y) can be written in 
the form 

a{x + 1, 7/) = a^i^x, y) a~^{x, y + l) a~^{x, y - I) a~^{x - 1, , (4.1) 

where we adopt the convention that 

a(x,0) = a(x,L2 + 1) = a(0, y) = a(Li + l,y) = 7, forallx,y. (4.2) 

The Eqns (4.1) can be recursively solved to express any operator a{x, y) as a product of 
powers of a(l,y). Therefore the group G can be generated by the L2 operators a(l,y). 
Denoting the rank of G (minimal number of generators) by g, this implies that 

g<L2. (4.3) 

In the special case of the linear chain, L2 = 1, we see that g — 1, and thus G is cyclic. 



From (4.1), a(Li + 1, y) can also be expressed as a product of powers of a(l, y), say 

a{Li + l,y)=Y[a{l,y'ryy' , (4.4) 

V' 

where the Uyyi are integers which depend on Li and L2, and which can be explicitly 
determined by solving the linear recurrence relation (4.1). The condition a{Li + 1, y) = 
I then leads to the closure relations 

a(l, y'Y'"'' = I, for alll < y < L2. (4.5) 

y'=l 

The Eqns (4.5) give a presentation of G, the structure of which can be determined from 
the normal form decomposition of the L2 x L2 matrix riyy'. This is considerably easier 
to handle than the normal form decomposition of the much larger matrix A needed for 
an arbitrary ASM. 

Even though this is a real computational improvement, the actual calculation, for 
arbitrary Li, of the rank of G is nontrivial. Even in the simplest case L2 = 2, it 
depends in a complicated way on the number-theoretic properties of Li. Details of 
this case are given in Appendix A. 

As to the square lattice where Li — L2 — L, using the above algorithm we find the 



following structures of G for the first values of L 

L = 2 : G = Z24X Zs, (4.6a) 

L = 3 : G = Z224 X Z112 X Z4, (4.66) 

L = 4 : G = Zqqoo X Z1320 X Zs X Zg, (4.6c) 

L = 5 : G = ^102960 X ■Z'l02960 X -^48 X Ziq X Z4. (4.6(i) 

This suggests, and it is in fact not difficult to prove, that for an L x L square 

g = L, for Li = L2 = L. (4.7) 



The idea of the proof of (4.7) is to use the fact that in this case, the matrix A has 
exactly L linearly independent eigenvectors of eigenvalue 4, of which all components 



can be chosen to be integers. We use them to construct L independent operators (none 
can be expressed as product of powers of the others) ^7fe,0<A;<L — 1, such that 

Ut = /. (4.8) 

This imphes that the number of generators is at least equal to L. Combined with the 
inequality (4.3), we get the result. 

To write the operators U]^ explicitly, we consider the eigenvectors nk{x,y) of A of 
eigenvalue 4 

A{x, y- x\ y')nk{x', y) = 4nfc(a;, y). (4.9) 

There are L independent solutions to Eq. (4.9), and a possible choice is to set 
{I < x,y < L, Q < k < L - I) 

nk{x,y) = {-iy-'[6{y-x-k)-6{y+x-k)-5{x + y-2L-2 + k)+6{y-x + k)]. (4.10) 

It is then easy to verify by using (4.9) and the relations (2.3) satisfied by the a{x,y), 
that the operators C//j defined by 

Uk = \[a{x,yr^^-^y\ (4.11) 

x,y 

all satisfy = I. Note that a strong form of independence of the eigenvectors nfc(x, y) 
must hold for the operators Uk to be multiplicatively independent: the eigenvectors 
must be linearly independent modulo the lattice {"^^i y, mx\y'A{x,y]x' ,y') : 
mx',j/' integers}, on account of the relations (2.3). The vectors (4.10) are easily seen to 
satisfy this condition. This completes the proof. 

When Li L2, the above construction is easily generahzed. If / = gcd[(Li + 
1), (L2 + 1)], it is easy to check that the matrix A has / — 1 independent eigenvectors 
of eigenvalue 4. Thus we obtain in this case that the rank QLixL^ satisfies the two 
inequalities 

f-l<gL,xL,<L2. (4.12) 

That the eigenvectors of integer eigenvalue play a particular role in the above proof 
should be clear. An arbitrary eigenvalue is generically an algebraic number, and the 



corresponding eigenvector has components which are also algebraic numbers, so that 
(4.11) becomes meaningless. However we show in the next two sections that one can 
construct toppling invariants from the spectrum of A, which yield the proper setting 
to generalize the above idea to any eigenvalue. The corresponding construction is 
applicable to any ASM which has a diagonalizable toppling matrix. 



5. A reminder of Galois theory. 

In this section, we recall the basic ideas of Galois theory, which can be used in 
any ASM, whatever its type. For illustrating the technique, we shall consider the 
familiar ASM defined on an L x L square lattice. The toppling matrix A is the discrete 
Laplacian, and its eigenvalues are easily determined and read 

^ ^ ^ 27rm 27m ^ , , 

Am,n = 4-2cos-^-2cos-^, l<m,n<L, (5.1) 

where we have defined 

A^ = 2(L+ 1) 
The order of the group G is given by 

L L 

|G| = n n (5.3) 

m=l n=l 

We define the transformations 



(5.2) 



(7-s(Am,n) = >^m',n', s coprime with TV, (5.4) 

where m',n' are integers, 1 < m' ,n' < L, such that cos(27Tm' /N) = cos{27rsm/N) 
and cos{2nn' / N) — cos{2Trsn/N). We see that these transformations act on the set of 
eigenvalues by permutation. Moreover, they form a group, and as 

(7s O as' = CTss' = CTs' ° CTs , (5.5) 

the group is commutative and isomorphic to the multiplication modulo N of the 
numbers coprime with N. Finally from the explicit expression of Xm,n, the actions 



of (Tg and o"_s arc the same. Thus the group of all ag is isomorphic to the group Z]^ 
of invertible integers modulo (those which are coprime with A^) quotiented by the 
subgroup {+1,-1} 

G&Il — {(Js ■ 1 < s < L and s coprime with N} ~ Z^/{±1}, (5.6) 

and is of order ^ip(N). We call the group (5.6) the Galois group of this ASM. 

It is instructive to group the set of eigenvalues Xm,n into orbits under the Galois 
group (5.6). As an example, let us take L = 3 or = 8. The Galois group consists 
of the two transformations s = 1 (the identity transformation) and s = 3. We find 
C3(Ai,i) = A3^3 and as{Xs^s) = so that {Ai,i, Aa^a} form an orbit under the Galois 
group. Gomputing the product of these two eigenvalues, we find Ai^iAa^a = 8. Likewise 
C3(A2,2) = A2,2 shows that A2,2 = 4 is an orbit on its own. (Note that because of the 
existence of some degeneracies in the eigenvalues, finding the orbit of Xm,n is not merely 
finding the orbit of (m, n) under a diagonal multiplication by all s. In the example at 
hand, to say that {Ai,i, As^s} is the orbit of Ai,i supposes that we have checked that 
the two eigenvalues are different complex numbers.) Doing the same calculation for 
the other eigenvalues, we can rearrange the product (5.3) giving the order of G as a 
product over the orbits (6 in this case) to find 

L = 3 : \G\ = [A11A33] [A22] [A12A32] [A21A23] [A13] [A31] 

(5.7) 

= 8-4 - 14 - 14-4-4- 100 352. 

We can do the same rearrangement for any value of L, writing \G\ as a product over 
orbits under the Galois group, 

1^1= n [ n Vn], (5.8) 

orbits O Xm,n^O 

and then each sub-product in square brackets is an integer. This follows from the fact 
that each square bracket is by construction invariant imder the whole of the Galois 
group, in which case the Galois theorem then states it must be a rational number. 
Because of the special form of the numbers Xm,n (they are algebraic integers), all 
square brackets must even be integers. Not only is each square bracket equal to an 



integer, but the Galois theorem asserts that none of them can be further spUt: the 
decomposition (5.8) is the maximal factorization of \G\ as a product of integers with 
the property that each integer is the product of a certain number of eigenvalues. 

Let us now put the above considerations in a more general perspective. A clear 
account of Galois theory can be found, for example, in [20]. 

Let P{x) be a finite degree polynomial with integer coefficients, and let us define 
the algebraic extension F[P] by adjoining to the field of rational numbers Q all the 
roots of P{x). (This kind of algebraic extensions are technically known as separable 
normal extensions.) Thus, F[P] is constructed in two steps: first we form the set of 
all rational linear combinations of the roots of P, and second, we promote this set to 
a field by adding whatever is needed to make it closed under multiplication. Hence 
F[P] contains all roots of P and all products of roots, but it is always possible to 
choose a finite number of base elements (finearly independent over Q) of which F[P] is 
the linear Q-span. The number of base elements is called the degree of the extension 
F[P], which has in general no relation with the degree of P. For instance the degree 
of F[x^ - 2] is 8, while that of ^[x^ - 1] is 2. 

Glearly, Q is a subfield of F[P], and it is thus a well-defined problem to look for 
automorphisms of F[P] which leave invariant every element oi Q C F[P]. The group 
of all such automorphisms is called the Galois group of F[P] with respect to Q, and 
noted Gal(-F[P]/Q). The effect of the Galois group on F[P] is to permute the roots of 
P, although not all permutations are allowed since they must be automorphisms. An 
equivalent definition of Gal{F[P]/Q) is the set of permutations of the roots of P{x) 
which preserve all algebraic relations among them. It can be shown that the order of 
the Galois group is equal to the degree of F[P]. We can now formulate the Galois 
Gorrespondence Theorem: there is a one-to-one correspondence between the subfields 
of F[P] and the subgroups of Gal{F[P]/Q). The bijection goes by associating a subfield 
M C F[P] with the maximal subgroup H C Gal{F[P]/Q) which leaves the elements 
of M invariant. Hence the bigger H, the smaller M. For instance M — Q is associated 
with H = Gal{F[P]/Q) by definition of the Galois group (hence an element of F[P] 
is in Q iff it is invariant under the whole Galois group), and M = F[P] is associated 
with the trivial subgroup H = {!}. 

Before closing this digression, we make a final comment about algebraic numbers 



and algebraic integers. All the elements of an algebraic extension of Q are algebraic 
numbers, which means that they all satisfy polynomial equations with integer 
coefficients. However, some of them may as well satisfy a polynomial equation with 
integer coefficients and with coefficient of the highest power equal to 1. These elements 
are distinguished in the extension and are called algebraic integers. For instance, in 
F[x^ — 2] = Q[-\/2], the number -\/2 is an algebraic integer whereas l/-\/2 is not. Also 
the eigenvalues of a matrix with integer entries are all algebraic integers. While the 
extension itself is a field, the subset of its algebraic integers is only a ring. In Q, this 
amounts to the usual distinction between the rationals and the integers. An immediate 
consequence is that an algebraic integer of F[P] which is invariant under the whole 
Galois group is actually an algebraic integer of Q, i-e. an integer of Z. 

Let us now see how these general ideas work if we take the extension F[det(x — A)] 
obtained by adjoining to Q the roots of the characteristic polynomial of the toppling 
matrix of an ASM. We choose the ASM defined in section 3 for an L x L lattice, with 
A the discrete Laplacian, but clearly the same analysis can be done in any model. 

We thus consider the algebraic extension Fl = F[det(a; — A^;,)] obtained by adding 
to Q the roots \m,n given in (5.1), 

_ , ^ 27rm ^ 27m ^ ^ ,^ ^, 

Am,n = 4 — 2 cos —^^^^ 2cos-^, l<m,n<L. (5.9) 

Since F^ contains all rational combinations of the Xm,n-, it contains cos for all 
1 < m < L since 

27rTO ^ 1 X /,- 

cos = 2 - — 2^ Xm,n- (5.10) 
n=l 

We also see from 

27rm 27TO 1 27r(m + n) 1 2t: im — n) 

cos -cos = -cos h - cos — , 5.11 

N N 2 N 2 N ' ^ ' 

that any product of Xm,n can be expressed as a rational combination of 1 and the 
numbers (5.10), which shows that any element of can be written as a rational 
combination of 1, cos cos . . ., cos ^j^. However, this writing is not unique as 



these L + 1 cosines are not independent over Q, but satisfy 

gco.^-(''^ + y + '^/'')=0. 0<a<^-l,0<6<p"-l, (5.12) 

for every prime power entering the prime decomposition of N . Using these relations, 
it can be shown that only ^ip{N) among the L + 1 cosines are independent over Q. 
Thus the extension Fj,, usually denoted by (5(cos^), has degree ^(p{N). A basis of 
Fl over Q can be taken to be {cos ^ : < m < lip{N) - 1}. 

To find the Galois group of Fl, we first note the following identities 

27rm 2tt 
cos-^ = T^(cos— ), (5.13) 

where is the m-th Chebyshev polynomial of the first kind. Equation (5.13) shows 
that the Galois transformation of cos^^, m > 2, is determined by that of cos^. 
The next step is to notice that, since the Galois group acts on the roots Xm,n by 
permutations, the Galois transformations of cos ^ must be in the set of cos^, 
1 < s < L. Hence suppose that o's(cos^) = cos is a Galois transformation for s 
between 1 and L. Being a group element, ag must have an inverse which we denote by 
as' = cr~^. It satisfies (Js'(cos^) = cos = cos^, hence ss' = ±1 mod N. This 
is impossible if s and N have a common factor, while if s is coprime with A^, there is 
a unique s' between 1 and L satisfying either s' — mod N or s' — —s~^ mod N. 
Thus the Galois group contains at least the transformations as for s between 1 and L 
and coprime with N. The number of such s is equal to ^(p{N), and since this must 
also be the order of the Galois group (equal to the degree of Fj,), one recovers the 
symmetry group of (5.6) 

GsIl ^{as : 1 < s < L and (s, N) = 1}. (5.14) 



We have completely identified the extension F^ and determined its Galois group. 
The last point concerns the ring of algebraic integers of Fl. This is a far more 
difficult question, and we just quote the result: the algebraic number w = qq + 
Sm=i cos of Fl is an algebraic integer if and only if all coefficients Qjtj 

are integers of Z. This ring of integers is usually denoted by Z(cos ^). In particular, 
as noted above, the eigenvalues Xm,n are all algebraic integers. 



Looking back at the decomposition (5.8) of \G\ as a product of orbits of eigenvalues 
under the Galois group, all the properties we mentioned easily follow from the above 
facts. The product of eigenvalues belonging to an orbit is by definition Galois invariant 
and thus must be an rational number. Since the eigenvalues are algebraic integers, this 
rational number must in fact be an integer. Finally the decomposition is maximal by 
definition of the orbits as the smallest subsets invariant under the Galois group. 

Similarly to the decomposition of G = Det A, we may consider the factorization of 
the characteristic polynomial of A in factors which are all invariant under the Galois 
group 

L 

Pl{x)^ l[{x-Xm,n)^ n [ n (^-Vn)]- (5-15) 
m,n=l orbits O Xm,n€0 

In this case, the Galois theory states that the polynomial within each square 
bracket, call it Po{x), has integer coefficients, and that the decomposition (5.15) is 
maximal, namely one cannot perform a further splitting into polynomials with rational 
coefficients (or in other words, all Po{x) are irreducible over the field of rationals). 

Examples of decomposition (5.15) for small values of L are 



P2{4-x) =a;2(x2 -4), 
P3(4 - x) = -x^ix"^ - 2)2(x2 - 8) , 

P4(4 - x) = x^(x - l)^(x + 1)2(^2 - 5f{x'^ -2x- 4)(j;2 + 2x - 4) , 

P7(4 -x) = -x\x^ - 2f{x^ - 8)(a;4 - 4^^ + 2f{x^ - %x^ + 8)^ 
(x^ - 16^2 + 32) (x^ - %x^ - 8a: - 2f[x^ - %x^ + 8x - if. (5.16d) 



(5.16a) 
(5.166) 
(5.16c) 



The number FiV) of factors in the decompositions (5.8) or (5.15) can be explicitely 
computed. The behaviour of F{V) is rather chaotic, as can be checked from Table I, 
where we give some of its values. It strongly depends on the prime decomposition of 
(L+ 1), and for large L, -F(L) increases linearly with L. An analytic formula for -F(L) 
is given in Appendix B. 



TABLE I: Number F{L) of irreducible 
factors in det{x — A). 



L 


F{L) 


L 




L 


F{L) 


L 


F{L) 


1 


1 


6 


16 


21 


74 


26 


112 


2 


4 


7 


19 


22 


64 


27 


112 


3 


6 


8 


28 


23 


122 


28 


84 


4 


12 


9 


34 


24 


88 


29 


192 


5 


18 


10 


28 


25 


90 


30 


88 



Finally, let us note that the Galois group (5.14) also acts on the L x L lattice 
on which the ASM is defined. This action is defined by as{x,y) = {x',y') where 
cos = cos and cos = cos The operator cTg permutes the lattice sites. 
In turn, this induces an action on the set of all configurations of the ASM by setting 

asC : as{zi) = z^^^i). (5.17) 

So we can also look at the Galois transformations either as automorphisms of the 
lattice or as automorphisms of the set of all configurations of the ASM. In fact, they 
even define automorphisms of TZ, because, although agC is not necessarily recurrent 
even if C is, there is a unique recurrent configuration equivalent to agC. 

6. Toppling invariants from eigenvectors. 

Using the basic results of Galois theory briefly recalled in section 5, we study here 
toppling invariants constructed from the (left) eigenvectors of A. This can be done for 
any ASM with diagonalizable toppling matrix, but for definiteness, we proceed with 
the ASM defined in section 4, for a square L x L lattice. 

Thus A is the discrete Laplacian, whose eigenvalues Xm,n were given in (5.1), and 
whose eigenvectors read 

sin^ sin^ 
vm,n{^: y) = 2^m " ^<m,n<L. (6.1) 

bill bill 

The normalization is purely conventional, and will be explained below. 



For each eigenvector Vm^i, we define an algebraic topphng invariant by taking its 
inner product with the z vector of an arbitrary stable configuration, and by evaluating 
the result modulo the eigenvalue Xm,n of Vm,n 

^m,n(C) = ^ Vm,n{x, y) Zx,y mod Xm,n- (6.2) 
x,y 

Since the eigenvalues and the eigenvector components are in general algebraic numbers 
in the algebraic extension F[det(a; — A)], equal to Q{cos ^) in this case, the congruence 
(6.2) must be understood in that extension. Moreover, modulo \m,n means modulo all 
multiples wXm,n with w in the ring of algebraic integers of Q{cos ^). 

Clearly under the toppling at site (xo.yo), we have Am,n — ^ ^m.n + 
Xm,nVni,nixo,yo), SO that A,„.„(C) is invariant provided Vm,n{xo,yo) is an algebraic 
integer for all XQ,yo. This was precisely the reason for choosing the normalization 
(6.1). On account of the fact that, for integer x, ^jj^ is a polynomial of 2 cos a with 
integer coefficients, the components of Vm,n are indeed algebraic integers of Q{cos ^) 
(see section 5). 

Thus the invariants Am,n are all valued in Z(cos ^) where the congruence (6.2) 
is to be taken. To give an explicit example, consider L = 3, = 8. The relevant 
extension is Q{cos^) = Q{\/2). Choosing for instance m = n = 1, one obtains 
Ai,i = 4 - 2V2, vi,i{x, y) = (1, V2, 1) ® (1, V2, 1), so that 

Al,l{C) = Zi^i + V2zi,2 + Zi,s + V2Z2,1 + 2Z2,2 + V2Z2,3 + ^3,1 

(6.3) 

+ V2z3,2 + Z3^3 mod 4 - 2V2. 

It can be seen that this invariant takes 8 different values by setting 1 = k for 
A; = 0, 1, . . . , 7, and all other Zx,y = 0. Indeed 8 = (4 + 2v^)(4 - 2V2) is the smallest 
positive integer equal to a multiple of 4 — 2y/2. Also recall from section 5 that Ai,i 
belongs to the orbit O — {Ai^i, As^s} under the Galois group, and that the associated 
irreducible polynomial is Pc'(x) = {x — Ai^i)(a; — As^s) = x^ — 8x + 8. So we have in 
this case that Po{0) = 8 is the number of values taken by the invariant Ai^i. 

It seems that we have defined invariants Am,n, one for each eigenvalue. However 
some of them are related by Galois transformations, and are thus not independent. 



Indeed if Xm.n and Xm^n' are in the same Galois orbit, the corresponding eigenvectors 
are related by a Galois transformation a, and so are the corresponding invariants, 
Am',n' = c(^m,ri)- So in fact we associate an invariant not to each eigenvalue, but 
rather to each orbit O under the Galois group, or equivalently to each irreducible 
factor Po{x) of the characteristic polynomial of A. We thus have F{L) (the number 
of orbits) algebraic invariants. Moreover, because of the normalization we chose in 
(6.1) for the eigenvectors, each invariant is linear in the height of the left top corner 
{x — y — 1) with coefficient 1: 

^m,n(C) = ^1,1 + . . . mod Xm,n- (6.4) 

Since 1 is not an eigenvalue, we see that none of the A^^n is trivial. However one 
cannot say in general the number of values they take. If one evaluates Am,n{C) at 
those configurations for which Zx,y = for all {x, y) ^ (1, 1) and zi^i = k ^ N, one sees 
that Am,n takes positive integer values, but because of the congruence, different /c's will 
give Am,n the same value. Denoting by kiaaxiin, n) the smallest positive integer equal to 
zero modulo Xm,n, Am,n will take exactly /cmax(?7i, n) different values. To determine that 
integer is not easy. Clearly an upper bound is given by the product of Xm,n with all its 
Galois conjugates, an integer by construction. This gives /cmax("^, n) < |P0(O)| where 
O is the orbit containing Xm,n- It can be shown that the equality holds whenever Xm,n 
is not the product of two integers belonging to two different subrings of Z{cos ^). In 
the other cases, the exact value of /cmax(^^5 n) depends on the factorization properties of 
Xm,n- [For instance, for L = 4, iV = 10, we have Q(cos |^) = Q{\/5). For Ai,i = 5 — \/5, 
we find /cmax = 10, while the associated value of Po(0) = (5 — a/5) (5 + v^) is equal 
to 20. The reason is that 5 — y/b — 2 ^~^ is the product of 2, an integer of Z and of 
an integer of Z{V5).] 

The next question is whether the invariants Am,n are linearly independent. This 
is easily answered negatively. There are F{L) > L invariants and we have shown that 
none of them is trivial, so that if they were independent, the group G would be the 
product of at least F{L) cyclic groups since it would have to contain the subgroup 
X m,n -^jfcmax (m,n)- (Remember from section 3 that toppling invariants provide additive 
representations of G.) This is of course impossible since, from section 4, we know that 
G is the product of exactly L cyclic factors. 



Even if the set of Am.n is not linearly independent, we can still consider the subset 
of independent ones, hopefully of cardinality equal to L, and see if those provide a 
complete labelling of the set of recurrent configurations. This too can be easily shown 
to fail. If they did, the total number of different values they take should be equal to 
Det A. But since Det A = Ylo \^o{^) I > the right number of values is reached only if we 
take all the invariants and if, for each of them, A;max("T., n) = |Po(0)|. This is clearly 
impossible from the linear dependence among some of them. 

Thus the invariants Am,n are neither independent nor complete. In view of the 
results of section 3, this is hardly surprising since the Smith normal form is not related 
to a spectral decomposition. A natural question is then whether one can find a complete 
set of algebraic invariants, analogous to the complete set of integer invariants /j defined 
in (3.11) in terms of the Smith normal form. It is indeed possible by using the following 
algorithm. 

We note that since the invariants involve congruences, there is no need to use the 
exact eigenvectors. We can define invariants by setting 



With no loss of generality, we can assume that the components of vx are in the ring 
of integers of some algebraic extension (otherwise A is simply rescaled), so that we 
can see the vectors vx as left eigenvectors of zero eigenvalue over a finite ring. We 
will call them modular eigenvectors (with zero eigenvalue). We look for a set of such 
vectors, as small as possible and with as large values of A as possible, in order to 
generate a complete set of invariants with the fewest invariants. To avoid irrelevant 
overall factors (if vx is a modular eigenvector modulo A, 2vx is a modular eigenvector 
modulo 2A), we may assume that at least one component of vx is coprime with A, 
which we can then choose equal to 1. The system (6.6) is overdetermined since it 
yields equations for only — 1 components of vx (one of them was set equal to 




(6.5) 



where we now take for vx a vector with the property that 



Vx ■ ^ = mod A. 



(6.6) 



1). The extra equation can be seen as a constraint on the values of A. By the Chinese 
Remainder theorem, it is sufficient to look for solutions of (6.6) for A a prime power. 
For if vx^ and vx^ are modular eigenvectors modulo Ai and A2 respectively, with Ai and 
A2 coprime, v — VX1X2 + vx^^i is a modular eigenvector modulo A1A2. Determining 
all independent solutions modulo prime powers and systematically using the Chinese 
Remainder theorem to obtain the smallest number of modular eigenvectors, we get a 
set of eigenvectors {vxi} where A^ is a multiple of Xi+i. ((6.6) implies that Ai, hence all 
Aj, divides the determinant of A.) The corresponding set of invariants {^Ia,} defines a 
complete set of toppling invariants. 

As an example, consider once more the 2x2 case. Let us choose the v^'s and the 
A's in Z, and assume 

^x{C) = zi + bz2 + czs + dz4 mod A, (6.7) 

(Zl Z2 \ 
j . The equation (6.6) implies 
Z3 Z4 J 

c — A — b mod A and d — Ab — 1 mod A, as well as 

86 - 16 = 166 - 8 = mod A. (6.8) 

For X — with p an odd prime, (6.8) is equivalent to 6 = 2 = mod p"*, which 
has no solution unless p^ = 3, in which case one finds one invariant, 

A3{C) ^ zi + 2z2 + 2z3 + Z4 mod 3. (6.9) 

For A = 2"^, the two equations (6.8) are trivially satisfied for 2"^ = 8, but cannot be 
satisfied for 2"^ > 16 (166 — 8 = —8 7^ mod 16). Thus the maximal value of A is 8, 
and one finds two independent invariants modulo 8 (6 = 1, —2) 

A8(C) = ^1 + ^2 + 3^3 + 3^4 mod 8. (6.10a) 
A'^{C) ^zi- 2z2 - 2z3 - Z4 mod 8, (6.106) 

Using the Chinese Remainder theorem, we obtain two invariants, modulo 24 and 8 
respectively: 



A24(C) = 8^3 (C) + 3A'g{C) = 11^1 + 10^2 + 10^3 + 5z4 mod 24, (6.11a) 



Agi^C) = zi + Z2 + 3^:3 + 82:4 mod 8. 



(6.116) 



As expected, one recovers the invariants (3.8) constructed out of the Smith normal 
form of A, 7i = 5^24 and I2 — SAg. Clearly, this is a general fact: choosing to work on 
Z, the above algorithm allows to compute the elementary divisors of A (the A^'s), as 
well as the relevant lines of the matrix needed to obtain the invariants (see (3.11)), 
since from (3.9) and (3.10), the rows of are precisely modular eigenvectors modulo 
the elementary divisors of A. However the algorithm is completely general regarding 
the ring of integers we want to work in. In particular, one can find a complete set of 
invariants with values in the ring of integers of any algebraic extension, something we 
will not pursue here since the algebraic invariants we would so obtain can be read off 
from the Z-valued ones (by decomposing the elementary divisors di & Z into a product 
of powers of prime ideals of the extension) . 

As a conclusion to this section, we found it rather attractive to use Galois theory 
in order to associate toppling invariants with irreducible factors of the characteristic 
polynomial of A, to the end of extracting information about G from the Galois 
decomposition of this polynomial. Although this method can be useful, the usefulness 
being dependent on the spectrum of A, it generally provides only partial information 
on G. 



7. Some related problems 

In this section, we comment briefly on two intriguing questions related to sandpile 
automata on finite square lattices. 

QUESTION 1. Wiesenfeld et al [15] have studied a deterministic version of the sandpile 
automaton considered here, in which the particle addition is always done at the central 
site of a (2L+ 1) x (2L + 1) square lattice. They observed, and it is easy to prove, that 
starting from an initial arbitrary configuration, once the transient configurations are 
gone, the system shows a cychc behavior, and the period of the cycle is independent 
of the initial configuration. The exact dependence of Tl on L is still quite puzzling. 



Wiesenfeld et a/ have numerically determined Tl = 4, 16, 104, 544, 146 248, 7 889 840 
for L = to 5, and found that Ti increases at a much slower rate than Det A, 
Tl ~ exp0.44L^ against Det A ~ exp4.67L^. The period for L up to 9 has been 
determined by Markosova and Markos [16]. It is easy to see that T^ is the order of 
the operator a(L+l,L + l)on the space TZ of recurrent configurations and so is the 
smallest positive integer such that 

a(L + l,L + l)^^ = 7. (7.1) 

If a{L + 1, L + 1)-^^ — I on TZ, we obtain that all the topphng invariants (3.3) must 
be equal to zero modulo 1 at the configuration which is zero everywhere except at 
the central site, where it is T^- Multiplying the invariants by Det A and defining the 
integer matrix E = (Det A)A~^, we obtain that T^ is the smallest positive integer 
satisfying 

EiXL+i,L+i) Tl^O mod (Det A), for all i. (7.2) 
Setting M — gcd {E^j^l^i l+i) ■ 1 < i < 2L + 1}, it is easy to see that 

rL = (DetA)/M. (7.3) 



Since Tl is independent of the initial configuration, we may choose a configuration 
which has the symmetry of the square. This symmetry is preserved under symmetrical 
topplings, and under addition of sand at the central site, hence we can restrict ourselves 
to T^sym, the subspace of recurrent configurations having the symmetry of a square. 
On a (2L+ 1) x (2L+ 1) lattice, a symmetrical configuration is completely specified by 
the heights in an octant having (-^+^K-^+^) sites. Instead of studying the symmetrical 
configurations on a square lattice, we can just as well study all recurrent configurations 
on an octant Ol with a new toppling matrix Agym given by 

{A^yu,)tj = , for i, j eOl, (7.4) 

f 

where the sum over j' is over sites which are related to j by the symmetries of the 
square (dihedral group of order 8). Even though it does not satisfy X]i(^sym)y > 



for all i, it can be shown that the general results of [6] hold for this case, so that for 
example 

where Gsym is the abelian group generated by the operators ai ,i & Ol, subjected to 
the relations (2.3) with A replaced by Agym- A simple calculation shows that 

DetAsym= n Vn, (7.6) 

m<n 

where the product over m and n is over odd values of m < n between 1 and (2L + 1), 
and with Am,n given in (5.1). For large L, |Gsym| varies as ~ exp0.58L^, so that 

Tf, increases at a substantially lower rate than |Gsym|- As in section 3, the structure 
of Gsym can be determined by computing the elementary divisors of Asym- We find for 
the first four values of L 



L = 


1 : 


■ Ggyra — 


Zl6 X Z2, 


(7.7a) 


L = 


2 : 


'■ Ggyra = 


Zi04 X Z4 X Z2, 


(7.76) 


L = 


3 : 


■ G^sym — 


^544 X Zs2 X Z2X Z2, 


(7.7c) 


L = 


4 : 


■ G^sym — 


-Z'l46248 X Zg X X Z2 X Z2. 


(7.7d) 



We see that the order of the largest cyclic group in Gsym is apparently nothing but the 
period Ti. In terms of the invariants (3.11) constructed from the Smith normal form 
of Agym, this would mean that the element of the first row oi A ^ corresponding to the 
central site is coprime with the largest elementary divisor d^^™, implying = d^^^. 
What is certainly true is that must divide d^^^ since a{L+l, L+1) generate a cyclic 
subgroup of Gsym, the largest one being of order d^^™. Although there are plausibility 
arguments that is actually equal to d^^^, the proof (or disproof) of this simple fact 
is still lacking. 

The question of the rank of Ggym can also be addressed. However here, an argument 
based on the eigenvectors of eigenvalue 4, similar to the one we used in section 4, does 
not work since the degeneracy of the eigenvalue 4 is only [^] + 1, with [x] the integral 
part of X. (It nevertheless implies that so many cyclic factors of Gsym have order 



multiple of 4.) But on the other hand, it is easy to see that there are exactly L + 1 
independent modular eigenvectors modulo 2 (of Agyui), which shows that the rank of 
Gsym is L + 1. (L + 1 is also an upper bound since all operators a{x, y) for y) E Ol 
can be expressed in terms of a(l, y) for y — 1,2, . . . , L + 1.) The independent modular 
eigenvectors have only one non-zero component, at any of the L + 1 positions on the 
principal diagonal of Ol: v^]^2(^iy) — ^x,k^y,k for 1 < A; < L + 1. This proves that 
each cyclic subgroup of Ggym has order multiple of 2. 

The fact that must be a divisor of |Gsym| already gives us a fairly good 
upperbound estimate of it. It is instructive to look at the values of IGsyml/Ti for 
small L. From the results of Markosova et al [10], we see that for L = up to 8, 
these values are 1, 2, 2^, 2^, 2^, 2^^ 2^, 2^2, 2^^ • 2701. Note that they show a rather 
irregular, non-monotonic behavior with L. They have a large power of 2 as divisor, 
but sometimes have other factors as well. 

The fluctuations in appear to be related to the appearance of degeneracies in the 
eigenvalue spectrum of Agym. For example, the fact that the eigenvalue 4 is ([^] + 1)- 
fold degenerate and the existence of L + 1 modular eigenvectors modulo 2 imply that 
Tl is a divisior of 2~(^+[^/2l^|Gsym|- Another accidental degeneracy in the spectrum of 
Gsym occurs when 2L + 2 is a multiple of 6. Writing 2L + 2 = M, we can easily verify 
that 



In the Galois factorization of |Gsym| this implies that there are some factors which are 
repeated. For example, for L = 8, the factor 2701 occurs twice in |Gsym|- While for a 
given value of L, the degeneracies, and the group structure of Gsym can be explicitly 
determined, it seems difficult at this stage to say much more about this "irregular" 
variation of Ti with L for general L. 

QUESTION 2. This concerns the structure of the identity configuration. Even the 
L X L square case shows nontrivial features. The identity configuration Cid is the 
unique recurrent configuration which is equivalent under toppling to the configuration 
with all heights zero. We have already noted that all toppling invariants Ii{Ci^) must 



^2£-m,2e+m — 




—m- 



(7.8) 



be zero. It also implies that C[d is the unique recurrent configuration such that 
Tj = ^j(A~^)jj^j(Cid) are integers for all i (by using the invariants (3.3)). This 
integer vector has the interesting interpretation that Ti is the number of topplings 
occurring at i when Cj^ is added to itself. 

The identity configuration shows complicated fractal structures. Some color-coded 
pictures of these may be found in [21]. We mention here two remarkable properties of 
the identity.* First, in the central area of a 2L x 2L lattice, there is a whole square 
where all sites have maximal height, equal to 3. The linear size of this central square 
was numerically measured for L up to 100, and goes hke 4L/5, so that the central 
square covers approximately 4/25 of the whole lattice. 

The second property, even more remarkable, is that the identity configuration on 
the (2L + 1) X (2L + 1) lattice appears to be related in a very simple way to that on the 
2L X 2L lattice. Indeed, if we divide the configuration C^^^^ into four equal squares 
and pull them apart by one lattice spacing so as to leave a cross in the middle, we get 
Cj^^^"*"^^ provided the heights on the cross are properly assigned. In obvious notation, 
if 



{B, areLxL blocks), (7.9) 



where the four blocks Bi are related by the symmetry transformations of the square 
(since Cj^ has that symmetry), then 




^(2L+l) 



/Bi Ri B2\ 

R2 ^mid ^3 

\B3 Ra bJ 



{Ri are 1 x L rows) 



(7.10) 



In addition the height at the center ^^mid is always 0, and the branches of the cross 
given by the Ri (also related by symmetry transformations) have a simple structure. 
The first instance of this phenomenon occurs when going from the 2 x 2 to the 3x3 



* We are indebted to Jean-Louis Ruelle for having run a program to compute the identity, and for 
having pointed out these two regularities. 



a 



(2) 
id 




C 



(3) 



id 



(7.11) 



lattices: 

/2 1 2\ 

10 1. 
\2 1 2/ 

More identity configurations are listed in Appendix C, for which this property may be 
checked (as well as the first one mentioned above). 

In addition, the array Tj defined in a previous paragraph also has the "cross" 
property (7.9-10). For instance, 
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2^ 



(7.12) 



Because the Cij and T arrays are related by = AT, the fact that both have the 
"cross" property fixes their values on the branches of the cross (for odd lattices). For 
instance, the rows in T corresponding to the Ri (see (7.10)) must be equal to the rows 
bordering them, as shown in (7.12) for 2L + 1 = 7. The first property mentioned for 
Cid) namely the existence of a big central square with constant values, does not hold 
for T. There is in fact in T a central square of constant values, but its hnear size 
cannot exceed 3 if Cj^ is to be recurrent at all (see (7.12)). As a final remark, we 
observed that the total number of topplings occurring when the identity is added to 
itself, i.e. equal to 160 and 235 for L = 6 and 7, grows like a power of L, 

with an exponent close to 4 (~ 3.9). 

Similar geometric structures are found when ASM are allowed to relax from special 
unstable states, say all heights equal to 4 [14]. These fractal structures are not well 
understood yet. 
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Appendix A 

In this appendix, we determine the structure of the group G for the ASM defined 
on the L X 2 lattice. In this case, the particle addition operators a{x, y) satisfy the 
relations 

a{x + 1, 1) = a^(x, 1) a-^{x - 1, 1) a"i(x, 2), {Al) 

and 

a{x + 1, 2) = a^{x, 2) a~^{x - 1, 2) a"^(a;, 1). (^2) 
Here 1 < x < L, and by convention we assume that 

a(0,7/) = a(L + l,|/) = /, for y = 1,2. (A3) 

Using the recursion relations (Al) and (A2) we can express a(x, 1) and a(x, 2) as 
product of powers of a = a(l, 1) and h = a(l, 2). Explicitely we have 

a{x,l) =a^'^ b-''^ (AA) 
a{x,2) ^a-"^ b'^'^, {A5) 

where /ix and Ux are positive integers satisfying the following recurrence relations 



{A6) 



with the initial conditions Ho = i^o = i^i = 0, and /i-i = 1. This can easily be solved by 
defining the new functions 

'mx = l^x- ^x, nj; = i^x + I'x- (AS) 
They satisfy uncoupled binary recurrence relations 

rrix+i = Srrix - rrix-i, {A9) 
rix+i = 5nx - rix-i, (AlO) 

with mo — riQ = 0, mi = ni = 1. The solutions read 
nix 
rix 

The first values of rux, Ux, Hx, are given below in Table II. 



TABLE II 



X 





1 


2 


3 


4 


5 


6 


7 


8 


nix 





1 


3 


8 


21 


55 


144 


377 


987 


Ux 





1 


5 


24 


115 


551 


2 640 


12 649 


60 605 







1 


4 


16 


68 


303 


1392 


6 513 


30 796 


I'x 








1 


8 


47 


248 


1248 


6136 


29 809 



One may note that nix are precisely the even terms in the standard Fibonacci 
sequence. 

To determine the structure of G in terms of the generators a and b, we note from 
{A3-5) that the only relations among them are 



1 

1 



2 

5 + V21' 



'3- Vs' 



5 - V21 



(^11) 
{A12) 



-VL+1 JjfJ-L+l _ J 



(AU) 



According to the discussions in sections 3 and 4, the structure of G is given in terms 
of the elementary divisors di, d2 of the matrix | ^^^^ ] fin(j ^^at 

Glx2 = X Zd^ , with di = ilJ^i+i - i^l+i)/d2 and d2 = gcd(/XL+i, z^l+i). (^15) 

As a function of L, the number d2 has an irregular behaviour, with very large sudden 
jumps. The only result we could firmly establish follows from the double inequality 
(4.12): ^2 > 1 if + 1 is divisible by 3, so in that case, G is not cychc. 

Rather surprisingly, the elementary divisors can equally be expressed in terms of 
mL_|_i and til+i by the neater formula 

di = ■mL+inL+i/d2, d2 = gcd(mi,+i, n^+i). (>116) 

Prom the relation (AS), Eqns (^15) and {A16) are clearly equivalent when rux and Ux 
are both odd, which is the case if x is not divisible by 3. That they are equivalent for 
X a multiple of 3 lies in the delicate fact that the largest power of 2 that divides rUx is 
the same as that dividing Ux- To see this, we consider the subsequences M(r) = msr 
and N{r) = nzr- We find from (All-12) that they satisfy the recursions 

M{r + 1) = 18 M(r) - M(r - 1), (^17) 
A^(r + 1) = llOA^(r) - A^(r- 1), (A18) 

subjected to the initial conditions M(0) = A^(0) = 0, M(l) = 8, and A^(l) = 24. Now 
the desired result that for each r, M(r) and N{r) have the same 2-potency (i.e. the 
largest powers of 2 dividing them are the same) follows from a more general result by 

the last author [22], which gives a formula for the 2-potency for any 'binary recursive 
sequence' d{r) defined by ^(0) = 0, ^(1) = 9, and e{r + l) = 2Re{r) + S e{r -1), where 
9 is an arbitrary integer, and R, S are odd integers. 



Appendix B 



In this appendix, we give the exphcit formula giving the number F{L) of irreducible 
factors (over Q) of the polynomial -Pl(x) for arbitrary L. In section 5, Pl{x) was the 
characteristic equation of the discrete Laplacian on a L x L lattice with open boundary 
conditions and was defined by (iV = 2(L + 1)) 

Pl(x)^ n U-4 + 2cos^ + 2cos^j. (Bl) 

■m,n=l ^ 

The roots \m,n = 4 — 2 cos — 2 cos ^fj^ of Pi belong to the extension Q{cos ^) and 
are permuted by its Galois group Gal/,, which acts on them by (see section 5) 

) = X,m,sn, for s e {±1}. iB2) 

As discussed in section 5, finding the irreducible factors of amounts to split the 
set of roots {A^^^} into orbits under Gal^, so that 

F{L) = \{Xm,n : l<m,n< L}/GalL|. {B3) 

By a classical theorem (see for instance [23]), the number of orbits of a set X under 
the action of a group G is equal to the average number of fixed-points of G in A, that 
is 



This allows to recast {B3) as 



2 ^ 



sez;,/{±i} 



From the expression of Xm,n, we have the obvious relations Xm,n = ^±m,±n = ^±n,±m 
where the signs are uncorrelated. The only subtle point of the analysis lies in the 
possible existence of other relations. For example, if L = 4 (i.e. A = 10), s takes the 



two values 1 and 3, we have that o"3(Ai,4) = A3, 2 is actually equal to Ai,4 = 4 even 
though (3,2) ^ (±1,±4) nor (±4, ±1) mod 10. This situation is however not generic, 
and can be proved to only happen for the pairs (m, n) with m + n = L + 1 (related to 
the eigenvalue 4 of the Laplacian). Thus if we leave them aside, we obtain the result 
that Xsm,sn — Xm,n if ^-nd Only if (sm, sn) — (±m, ±n) or (±n, ±m) mod N. Each of 
the L eigenvalues Xrn,L+i-m — 4 clearly forms an orbit on its own under the Galois 
group and we obtain 

F(L) = L+— — \{1 <m,n < L and m + n^ L + 1 : 

K€Z*^/{±1} 

{sm, sn) — ±(m, ±n) or ± {n, ±m) mod N}\. 

{B6) 

The actual calculation of (-B6) being straightforward but somewhat lengthy, we 
only quote the final result. Let N = 2(L + 1) have the prime factorization given by 

p 

Then the number of irreducible polynomial factors of Pl{x) equals 

F(L)^- + 6 + -n ^ _^ + 2(2a2-l)n(2ap + l) 

P P^2 

+ H (2ap + l)-(7a2 + 5) J](ap + 1) (^g) 

p=l mod 4 P7^2 

+ [Y[i2ap + l)-l]-5{a2 = l). 

It can be shown that the average growth of F{L) is linear in L. Some particular values 
are given in Table I, section 5. 



Appendix C 

We hst here the identity configurations C^^\ discussed in sections 3 and 7, for 
L — A and 5 and L — 10 and 11. 



/2 3 3 2\ 

3 2 2 3 

3 2 2 3 

V2 3 3 2/ 



(10) 
id 
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3 2 12 3 

2 10 12 

3 2 12 3 
V2 3 2 3 2/ 
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Appendix D 



In [17], Lee and Tzeng have considered the following one-dimensional deterministic 



sandpile model. The L sites of a linear chain, labelled to L — 1, are assigned a height 
variable Zi taking positive integer values. A site j becomes unstable if the difference 
Zj — Zj-\-i exceeds some integer N, considered as a parameter. In this case, the site j 
looses N grains of sand which fall on its right nearest neighbours. They take open 
boundary conditions Zi — for all i > L, so that, under the toppling at site j, the 
sand is redistributed according to Zi ^ Zi — Ay , where the toppling matrix is the lower 
triangular matrix 



The deterministic dynamics is defined by dropping, at each time step, one grain of sand 
on the leftmost site i = and by letting the system relax. The number of recurrent 

configurations is equal to N^, also equal to the determinant of A. Note that this model 
is a non-abelian sandpile model as the toppling condition is not local. 

It was shown in [17] that for = 2 and 3, the recurrent configurations can be 
labelled by the values of a toppling invariant, an integer function taken modulo A"^, 
(see also [24]). They wrote down the invariant when A" = 2 for all values of L, as well 
as when A" = 3 for some particular values of L. We show here that our general method 
allows to prove that one invariant is sufficient to label all recurrent configurations, for 
all A^, and we will compute this invariant for A^ = 2, 3, all values of L. 

That only one invariant gives a complete labelling of the recurrent configurations 
amounts to prove that the elementary divisors of A are di — and di — 1 ior i > 2. 
Equivalently, it amounts to show that among all the invariants (3.3) constructed from 
A"-"^, only one is independent. The matrix A"-*^ is easily obtained and read 





L-l 



where the entries satisfy the recurrence relations 



NCk = Ck-l + Ck-2 + • • • + Cfc-iV, 

Co = — , c_fc = for all k > 1. 



{D3) 



We thus obtain the L invariants 





k=o 



The first invariants read Qq = j^zq mod 1, Qi = -^zq + j^zi mod 1, Q2 = ^^t-zq + 
772-^1 + jfZ2 mod 1, so that we have 



NQi = Qo mod 1, 

NQi ^(N + l)Qo mod 1, NQ2 ^ (N + l)Qi mod 1, 



{D6a) 
{D6b) 



showing that Qo can be expressed in terms of Qi, itself expressable in terms of Q2- It 

is not difficult to show recursively that this pattern continues to hold, namely there 
exists, for all i between 1 and L — 1, an integer Xi coprime with N such that 



NQj — XiQj-i mod 1, for all I < j <i. 



(Prom (D6) we have Xi — 1 and X2 — N+1.) Equation (D7) implies that all invariants 
Qq,Qi, ■ ■ ■ , Ql-2 can be expressed in terms of Ql-i which is thus the only independent 
one and which provides a complete labelling of the set of recurrent configurations. 

To compute the invariant Ql-i, we note that the coefficients satisfy in fact a 
recurrence relation of order — 1 : 

Nck + {N- l)cfc_i + . . . + Ck+i-N = NcN-i + iN- l)cN-2 + . . . + CO = 1. {D8) 
For = 2 and 3, the solutions to (D8) read {k > 0) 



^ = 2 : Cfe = i[l-(^)^+i] 



" 6 12 



{D9a) 
{D9b) 



Using these values in Ql-i of (D4) yields the desired result. For higher values of N, 
the coefficients are given by the generating function 



fN{t)=j:ci''¥={i-t). 



^N+l 



{N + l)t + N 



-1 



{DIO) 



k=0 
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